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Abstract. We provide explicit formulae for highest-weight to highest- weight correlation 
functions of perfect vertex operators of U q (sl(2)) in arbitrary integer level I. They are given 
in terms of certain Macdonald polynomials. We apply this construction to the computation 
of the ground state of higher spin vertex models, spin chains (spin 1/2 XXZ) or loop models 
in the root of unity case q = —e^™^ e+2 \ 



1. Introduction 

In the seminal paper [13], Frenkel and Reshetikhin showed that correlation functions of 
g-deformed vertex operators (VOs) associated to quantized affine algebras satisfy a set of 
holonomic g-difference equations which are a g-deformation of the Knizhnik-Zamolodchikov 
equation [20] (gKZ). The analysis is performed for generic value of the level I of the affine 
algebra, but also works for positive integer level provided the obvious modifications are 
made. The most important one is that the paths in the Weyl chamber describing the various 
conformal blocks (correlation functions of vertex operators) have to be restricted to the 

"Weyl alcove". In all that follows we focus on the simplest algebra, U q (si(2)). In this 
case, the integrable irreducible highest weight modules are characterized by a spin s G |Z + 
(defined e.g. as the spin of the U q {s\(2)) representation of the top degree part, where U q (sl(2)) 

is the horizontal subalgebra of U q (si(2))), which in level I only exist when s < £/2. Each 
vertex operator corresponds to a step of the path and is itself characterized by a spin j £ 
{1/2, ... , £/2}. The present paper is entirely dedicated to the case of so-called perfect vertex 
operators, that is, j has its maximal value j = £/2, for which there is exactly one possible 
step, namely, s — » 1/2 — s. This implies that various simplifications occur, and we expect this 
unique conformal block to be particularly simple; and indeed we provide an explicit formula 
for them. 

The reason to revisit this somewhat old subject is the observation made in [7] that cer- 
tain solutions of gKZ also provide the ground state entries of integrable models at special 
values of the deformation parameter q. In fact, this idea can be traced back to [25], where 
Reshetikhin obtained solutions of gKZ using an "off-shell Bethe Ansatz" and observed that 
Bethe equations appeared in the semi-classical approximation, and then pursued in [26], 
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where the Knizhnik-Zamolodchikov equation in the quasi-classical limit was found to pro- 
duce eigenvectors of the Gaudin model. Here, we are only concerned with one especially 
simple solution of qKZ - with appropriate normalization, a polynomial solution, which pro- 
duces one eigenvector of the higher spin U q (sl(2)) (inhomogeneous, twisted) transfer matrix. 

Note that in earlier work [7] (see also [6, 24]), only level 1 solutions were used, whose 
explicit form was already known and can be found e.g. in [17]. In the present paper, we 
obtain a novel expression for the correlation functions of level i perfect VOs, and show that 
once specialized, they provide the ground state entries of the integrable spin £/2 chain when 
q = —e~ l7T /( i+2 \ We also establish the connection with the loop model of [c I] and thus prove 
Conj. 1 in it, which is concerned with the degree of the ground state entries as polynomials. 

Also note that a different procedure to produce (arbitrary level) solutions of the qKZ 
equation is to compute finite temperature correlation functions of the type tr(x d $ . . . $), 
see e.g. [17]. However the resulting integral formulae are of a more complicated nature than 
those we consider here. In particular, they only make sense for \q\ < 1 and it is not clear how 
to continue them to \q\ = 1. Here we restrict ourselves to the zero temperature correlation 
functions of the type (0| $ . . . $ |0), though tr(x d $ . . . $) should be computable along similar 
lines. 

The plan of the paper is as follows: in section 2 we review the construction of U q {sl{2)) 
currents in terms of bosons and parafermions and then proceed to build perfect vertex oper- 
ators. Section 3, the core of the paper, contains the derivation of the correlation function of 
these vertex operators. Section 4 then reinterprets this correlation function as the eigenvector 
of an integrable transfer matrix in various contexts (spin chain, loop model, supersymmetric 
lattice fermions). Finally, various appendices provide additional technical details. 



2. Construction of U q {sl{2)) perfect vertex operators 

2.1. The algebra U q (sl(2)). Our reference for the quantized affine algebra U = U q (sl(2)) 
is the book [17] as well as the paper [16], whose conventions we follow. The Chevalley 
generators are denoted by E^F^Ki, i = 0, 1, to which one must add the grading operator d: 
we choose the homogeneous gradation, i.e., E (resp. F ) has degree +1 (resp. —1) and all 
other generators have degree 0. In particular there is a horizontal subalgebra U\<zU which 
is generated by E\,F\,K\. 

We also use in what follows Drinfeld's realization of U in terms of currents [9]. 

Until specified otherwise we are in the regime \q\ < 1. 

We shall consider the spin £/2 evaluation representation p z acting on V z = C e+1 with 
standard basis v b , b — 0, . . . , I. It is defined by: 

p z {Fx)v h = [6]u6_i p z {F )v b =[£- b]z- 1 q" +2 v b+l 

(2.1) P*(£iH = [t- b]v b+1 p z {E )v h = [bjzq-^Vb-x 

PziKjv,, = q 2b - £ v b Pz (K )v b = q e - 2b v b 

where [n] := (q n — q~ n )/{q — q^ 1 )- In order to make d act, one must consider z a formal 
variable, in which case d = z-^ + A, where A is a constant to be fixed below. When restricted 
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to Ei,Fi,Ki, it is an ordinary irreducible representation of Ui = U q (sl(2)) of spin £/2, and 
we simply denote the underlying space V. 

We shall also need the R- matrix R(z\/ z 2 ) acting on a pair of such representations V Zl ®V Z2 . 
The bar denotes the fact that we choose for now a particular normalization of the i?-matrix 
where it is the identity on the highest weight vectors. Other normalizations will be considered 
below. Equivalently, consider R(z) = VR(z) where V permutes factors of the tensor product. 
We then have the following formula: 

(«) *w = E LT 

j=0 r=j+l 

where Pj is the projector onto the spin j irreducible subrepresentation of V <8> V w.r.t. the 
C/i-action. 

Finally, we shall denote collectively by He any level i representation space of U, with 
a : U ->• C(H t ). 



2.2. Level t U q (sl(2)) currents. The results which are summarized here are based on 
the work of Ding and Feigin [8] (see also [4]). For the sake of completeness and because 
some details do not appear in this article we give a detailed introduction to the subject in 
appendix A. 

Consider the bosonic operators: 

[2n] [in] 

(2.3) [a n , a m ] = 5 m _„ [a n ,(3} = 5 nfi n,meZ 



n 



with grading 



[d,a n }=na n [d,/3] = -^ 



They act on the direct sum of copies of the bosonic Fock space %f = © ftgZ (lIm>o a -m \^)b) 
which differ by the eigenvalue of do: a \h) B = h\h) B . (3 only appears via its exponential, 
which is such that e ±/3 \h) B — \h ± 1) B . The grading is given by d\h) B = —^h 2 \ h) B . When 
I — 1, this is a realization of U, and therefore it is denoted by T-L\. 

We build the following current operators: 

e (z) = ^ + (z)e 2 ^z ao/i e En> ° 9 ^w^V E " >0<? ^S 2 "" 

(2-4) ^) = ^-^ e -2/9 z -a /, e -E„>o9^^f ^gE^o^S 2 "" 

fe ± (^) = g ±a e ± (9-9" 1 )En a ±n 2:=Fn 



4 T. FONSECA AND P. ZINN- JUSTIN 

where the £ (z) are parafermions, see appendix A. A full list of relations they satisfy is given 
in (A. 2), as well as the relation with Chevalley generators. Here we list a few: 

k+{z)f{w) = Z ' q2+ ] W f(w)k + (z) 
q 2 z — qzw 

^2 5) e(z)e(w) [z — q 2 w) = e(w)e(z) (q 2 z — w) 

f(z)f(w) (q 2 z - w) = f(w)f(z) (z - q 2 w) 

(q - q- l )zw[e{z), f(w)} = 6 (V^) k + {q^z) - 5 (g^) AT(^) 
where 8{z) = E ie z^- 

These operators act on the tensor product of £ copies of the Fock space Hi, wich is a level 
£ representation of U. A state |A;i) B <S> ■ ■ • <S> \kg) B is denoted by \ki, . . . , kg). As we shall see 
in Section 3, when ^ = or 1 for all 2, \ki, . . . , k() is a highest weight vector, with highest 
weight given by the sum k = ^ k%\ for fixed k = ki, these states are interchangeable and 
we denote any of them by \k). 

Remark: for an alternative construction of these vertex operators, see [5]. 
2.3. Intertwining relations. Let be a type I vertex operator, i.e., an intertwiner 

whose normalization will be fixed later. 
The intertwining condition writes: 

(2.6) (a ® p g )(A(x))$(z) = §{z)a{x) 
for all x £ U, or explicitly for the Chevalley generators: 

(a(Ei) ® 1 + a(Ki) g> p t (Ei))$(z) = <$>(z)a(Ei) 

(2.7) (a(Fi) <g> Pz {K7 l ) + 1® Ps (FJ)$(z) = $(z)<r(Fj 

(a(Ki) ® p s (Ki))$(z) = $(z)*(Ki) 
plus the intertwining condition for d which will be discussed separately below. 

In what follows application of a will be implicit. Using the form of p z , cf (2.1), we can 
expand in components: $>(z) = ^2 b=0 $> b (z)v b with 

§ b (z)E - E $ b (z) = [b + l}zq- e - 2 K <S> b+1 (z) 
& b (z)Fi - q e - 2b F^ b (z) = [b+ l]$ b+ i(z) 

(2.8) Mz)Ei - Ei$ b (z) = [£-b+ qK&t-.tiz) 
$ b (z)F - q 2b - l F^ b {z) = [l-b+ l]z-V+ 2 $ 6 _ 1 ( Z ) 

q 2b - £ K^ b (z) = <Z> b K! 

Since He is a level £ representation, we have KqKi = q 1 . 

We are mostly interested in operators F\ and Eq; expanding the lowering current f(z) in 
modes: f(z) = J2nez fn z ~ n ~ 1 , we have f = Fi, f\ = E Ki, so that we rewrite the two first 
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equations in (2.8): 

**(s)/o - g'- 26 /o$ 6 (^) = [b + l]$ b+ i(z) 

' Mz)fi-q 2b - e fiMz)=q 2b ~ e [b + i]z® b+ i(z) 

This means starting from $o(z), one can build the by iterated g-commutator. We obtain: 



(2.10) k=0 J 

/ i \fe„(6-fe)(£-6+l) 

"* ^ [*]![*-*]! /l<W/l 

where [n]! = [n][n — 1] . . . [1]. 

Finally, we write / = § f(w)dw, fi — § f(w)wdw, so that: 
(2.11) 

/C b fin,, b C_1 U-Js(€-6+l) 
■ ■ • / II E [fc ]![?-fc]! • • • /H)^ WH + x) • • • /(«*) 



i=l fc=0 
6 , /. 



/r " in. rlnn f— 1 \k n (b-k)(l-b+l) 

■ ■ ■ j II ^ E ' /W • • • /K)*oW/K + • • • /K) 

The summation over will be performed in section 2.4.2. 

2.4. Explicit formulae for the perfect vertex operators. So far we have not used the 
fact that twice the spin of the VO is equal to the level. For so-called perfect VOs, that is for 
rank 1 case precisely the VOs of maximal spin £/2, we also have the following form of the 
highest weight entry in terms of the bosonic field: 

(2.12) «&„(*) = e ^(- z )°°/* e E^^V^ ¥ W z - n 

In order to check the consistency of this Ansatz, we shall require certain identities. 

2.4.1. Commutation relations. In order to compute the commutation relations between the 
operators e(w), f(w) and &o(z), it is convenient to define the normal ordering :...:, which 
pushes the negative modes to the left, that is: 

:a. m a n : \ :a f3: \ _„ 

:a n a_ m . j ~ a -°» :/3a : J " ^° 

where n and m are positive integers. The normal ordering has the following important 
properties: 

Lemma 2.1. Let O be a product of N operators O = 0\ . . . On- If OiOj = gij : OiOj : , 
then O II,. , .</,,, - O: . 

Lemma 2.2. Let O be an operator such that :0: = e s3 /(a )e^»>» d » a -"e^»>« c " a " , where 
g, c n and d n are scalar constants. Then 

B (k\ :0: \k) B = f(k) 
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if g = and it vanishes otherwise. 

In the more general context where :0: = <^) i e 9ia fi(ao)e^ n >° di ^ a -™e^ n >° Ci >™ an this property 
is conserved: 

(kx,...,ki\ :0: \ki,...,k e ) = 

i 

if gi = for all i and it vanishes otherwise. 

For convenience, we split the currents in bosonic and parafermionic components, that is, 
e(w) = £ + (u>)e(u>) and f(w) = £~(w)f(w). 

Proposition 2.3. The relation between the products and the normal ordered products is: 

<? 2 f \q 

%( Zl )<Z> (z 2 ) = {-z x )i v ^-22-: $0(^1)^0(^2): 

f(«,)$ (z) = _J- r :f(«,)$ (z): 

$ (^)f( w ) = _J_:f( w )$ (^): 

g £ w — 2; 

c(w)$ (-2) = (w - z) : e(w)$ (-2) : 
® (z)t(w) = (w — z): t(w)%(z) : 

Remark 2.4. In t/ie previous computations, we perform infinite sums which are convergent 
under the following conditions: 

\q 2 Z2\ < \z\\ \q z\ < \w\ < \z\ \z\ < \w\ \w\ < \z\ 

for each expression in proposition 2.3, respectively. 

From these expressions we find immediately: 

(2.13) aw — z 

As proven in appendix A, the parafermions do not interact with the bosons. Therefore, in 
these formulae we can replace f(w) with f(w). 

From the first identity we conclude that [e&, $0(2)] — 0. For k = 0,1, this coincides 
with the intertwining equations in (2.8) for E\ and F at b = 0. We of course also have 
Ki§q = q e § Ki, using K\ = q a °. Finally, we have: 

(2.14) x d %{z)x~ d = x- E/ % Q {x- l z) 

which coincides with the intertwining condition for d on condition that one set A = £/4. 1 



1 Note that in Conformal Field Theory (corresponding to q — > 1), this would be interpreted as a con- 
sequence of the conformal dimension £/A of $>o(z). The q-deformation breaks the conformal symmetry by 
marking the two points and 00, leaving only the residual scaling symmetry z xz. 
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2.4.2. Closed expression for the $ b (w). Equation (2.10) provides us with two expressions for 
$5(2). Performing the summation results in: 

Lemma 2.5. The following expression holds: 

Widwi $oO)/(«t) • • • f(wb) 



<$> b (z) = a b j> ... J>Y[ 
J i=i 



2ni Hi(wi - q e z) 



where <x b is a constant given by: 



[&]! v ' ^ * ' ^ [b]\[e-b]\ 

Proof. We shall use the first expression of &b{z) in (2.11); the second expression would lead 
to the same result. We have a sum over several terms. We commute all / to the right, so 
that each term becomes: 

/M • • • f(w k )$ (z)f(w k+1 ) . . . f(w b ) = ^- ki f W3 _~ t Z U o(z)fM . . . f(w b ) 

LLj<k\ w j 9 z ) 
FT ^u(Q Wi — Z) FT ..,(«); — n i z) 

- *° <M ■ ■ ■ ' « 

Our purpose is to compute the expectation value of an operator, which contains the 
operator Q (z)f(wi) . . . f{wb)- By lemmas 2.1 and 2.2, this will make appear one term per 
each pair of operators. From §Q(z)f(wi) we get a pole g i J._ z ■ The contribution of the pairs 
f(wi)f(wj) is more complicated and we need to consider the decomposition of parafermions 
in modes f k (wi), which we explain in appendix A. Here we only need a couple of properties. 

The relations (A. 2), which define U q (sl(2)), imply that 

±± Wi — Wj 

Kj J 

is symmetric in the exchange of the Wi. This expression decomposes as a sum of terms of 
the form f n (w\) . . . f lh {wb) Yli<j(l 2w i ~ w j)i which each satisfy by equation (A. 23) 

,f >x) . . . f"(w b ) l[(q 2 w t - Wj ) = P( Wl , ...,w b ): f^(w 1 ) . . . f>(w b ) : 

i<j 

where P is some multivariate polynomial. We conclude from the definition of normal ordering 

2 

that f(wi) . . . f(w b ) Yli<j is a symmetric function of the {w±, . . . , w b } without poles 

of the form {wi — (wj)^ 1 . 

In summary, if we write f b (z; Wi, . . . , w b ) = Ylj<i(i ew j ~ z ) Ylj>i( w j — ( f z ) (which is linear 
in each variable Wj), then the quantity that we want to compute is: 

1 — 1 k — ^"^J 
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where S(z; Wi, . . . , w b ) is some symmetric function of the {w±, . . . , w b } without poles of the 
form (wi — (Wj)^ 1 . 

We now want to prove that 



b i -t \ k h(V—h-i--\\ b 



^ [fc]![6-fc]! fb(^w 1 ,...,w b )=a b Y[w i 

up to some terms which vanish when integrated, i.e., which are divisible by (q 2 Wi — iUj+i). 
They produce a zero contribution because the integral becomes skew-symmetric. 

This is true when 6 = 1: 

f?(z; wi) - q e fl(z; Wl ) = (wi - <fz) - q\q l w x - z) = (1 - q 2l ) Wl 

Suppose that this is true for 6 — 1, then: 

fc — /c — y ^k j^>k 

_]\k-iq(k-i)(e-b+2) 



k — 1 j "^k j^k 

-^p(™& - q z) [[Wi - (q w 1 - z) [[wi 



by the induction hypothesis. The last term can be rewritten: 

i=l \j<i J \j>i+l 

which vanishes when integrated. 

The constant a b is given by this recursion. □ 



2.5. The R matrix. Consider next the operator R(z-i/ z 2 )Q>{zi)$(z 2 ) : He — > 7ie.®V Z2 ®V zl , 
with the normalization of the i?-matrix to be adjusted below. R{z,\j 'z 2 )Q(zi)Q(z 2 ) is an 
intertwiner by the defining properties of R and $. We now use the fact that $ is a perfect 
VO, i.e., that it is of (maximal) spin £/2 where t is the level. This means that each highest 
weight irreducible representation (say of spin s) in He, when tensored with V z , is sent onto 
another unique irreducible representation (namely, of spin £/2 — s), and the same when 
tensored twice. Such an intertwiner is therefore unique; so that it should be proportional to 
$(2:2) $(2:1). In order to fix the proportionality constant, we compute 



r( Zl /z 2 )& \ Zl )& \z 2 ) = &°\z 2 )&°\ Zl ) 
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with 

rlz) = _.-</, («■«•■ f). 

1 ' (<i 2 A; 9 4 )oo (9 2f+2 ^l 9 4 )oo 

and now choose (as in [16]) 

R(z) = r(z)R(z) 

where R(z) is the identity on the tensor product of highest weight vectors, in such a way 
that 

(2.15) R(zi/z 2 )$(z 1 )<l>(z 2 ) = $(* 2 )$(*i) 

3. Correlation functions of U q (sl(2)) perfect vertex operators 
In this section we explain how to compute the correlation function in even size L = 2n: 

^....fcJZl. • • • > *2n) = (k\ $ 6l (^l) . . . <£v(^2n) 

\&6i 6 (^i, . . . , z 2n ) is nonzero only when the neutrality condition X]?=i &t = ^ satisfied. 
The final result will be a multiple contour integral. 

3.1. Integral formulae for correlation functions. By lemma 2.5, the product of vertex 
operators $^(21) . . . &b 2n ( z 2n) can be expressed as a multiple integral containing 2n opera- 
tors &o(zi) and in operators f(wj). Then we can use lemmas 2.1 and 2.2 to compute the 
correlation function. 

3.1.1. Correlation functions of q-parafermions. Take a product of in currents, multiplied as 
before by some appropriate prefactors, and do a mode decomposition: 

(3.1) T = f(w 1 )...f(w en )H^^= y: fw-f»wn— 

■ LX Wj — Wi xx Wi — Wj 

Kj J »1 ,»•>*/« Kj J 

Equation (A. 23) implies that the correlation function (0| q 2fn ^J r |0) is a symmetric polynomial 
in the u>j. 2 Moreover, as a consequence of theorem (A. 2) for its parafermionic part, it satisfies 
a wheel condition. 

This wheel condition can be defined in a general setting as follows. Let tM, Qm be two 
scalars such that q 7 ^ 1 ^ 1 = 1. A symmetric polynomial is said to satisfy the (r, fc)-wheel 
condition, if it vanishes whenever we set the first k + 1 variables such that ^± = t^q M for 
all i < k and Yli s i ^ r — 1. Here we find that (0| q 2tn ^T |0) satisfies the wheel condition 
when we set r = 2, k = i and tM = q~ 2 - 

For each term on the decomposition (3.1) we can compute the effect of normal ordering: 

r >i) . . . rw = i'" . . ■ , wtn) : r >i) . . . r-K) : 

Now, we ignore everything but the terms in e a . On the one hand, we have e 2£n ^ = . e na . 
On the other hand, each parafermionic mode f k is proportional to 1 ® e ~ a &)j>k 1- 

Then, by lemma 2.2, each mode k should appear exacly n times. 



'It can be shown that the coefficients of the polynomial are Laurent polynomials in q. 
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Each time we have a repeated mode (that is, q = 6j) we get a term of degree 2, more 
precisely q 2£ ~ 2ei (wi — q 2 Wj)(q 2 Wi — Wj). Otherwise, we get a term of degree zero. Therefore, 
the degree of (0| e ^J- |0) is £n{n — 1). The same analysis can be used to prove that 
(EL w i) (k\ e 2in P \k) is a symmetric polynomial of degree £n(n — 1) + n(£ — k). 

We can do a more detailed analysis and obtain the dominant monomial. Let Y^ n be the 
staircase Young diagram, corresponding to n steps of £ x 2: Y^ n = (2(n — 1), . . . ,2(n — 
1), 2{n - 2), . . . , 2(n - 2), 2, . . . , 2, 0, . . . , 0). For example: 



Y, 



3.3 



2.1 



(fc) 



More generally we define a modified staircase Young diagram Y £ n 
Ye n by adding £ — boxes vertically at each step. For example: 



Which is obtained from 



Y 



(0) 



3.3 



(i) 



3,3 



Y 



(2) 



3.3 



Yo 



(3) 



3,3 



where the doted boxes correspond to the boxes added to Y^ n . 

Recall that a monomial corresponding to some partition A = {Ai, 



A at} is defined by 



w 



Proposition 3.1. The leading monomial of (Y[i w i) (k\ e 2ln P F \k) is: 



In the case k = 0, the expression simplifies: 

(0| e 2 ^ |0) = /«("-i)+n(a) ([£}\) n w Y ^ + lower terms. 

Sketch of proof: multiply (Y[ i Wi) (k\ e^^T \k) by the Vandermonde determinant, and try 
to maximize the degree of w±, then try to maximize the degree of W2 and so on. 

Call a Young diagram (r, fc)-admissible if Aj — Xi+k > r for all i. Then, it is clear that 
a staircase (r, /c)-admissible Young diagram has the minimum number of boxes possible. If 
we set r = 2, k = £ and we restrict ourselves to diagrams of lenght 2n the only diagram 
with £n(n — 1) boxes is exacly Y^ n . For the case of modified staircases, the situation is 
different: there are several (2, £)-admissible Young diagrams with £n(n — 1) + (i — k)n boxes, 

(k) 

the smallest of them being Y £ n . 

Macdonald polynomials [ ] form a basis for the symmetric polynomials depending on two 
extra variables which we call qu and tu- The only result related to Macdonald polynomials 
that we need for our purposes is the following theorem: 

Theorem 3.2 (Feigin, Jimbo, Miwa and Mukhin [11]). Let qu, tu be such that q^ 1 ^ 1 = 1. 
Let V be the vector space of symmetric polyynomials on N variables satisfying the (r, k) -wheel 
condition. And let M. be the vector space spanned by the Macdonald polynomials given by 
(r, k) -admissible Young diagrams and the Macdonald parameters qM, tu- 
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Then V = M. 

As a direct consequence we can compute the correlation function: 
Theorem 3.3. 

(0| e 2 ^f( Wl ) . . . f(w in ) |0> TT Wj - q2Wl = 7 ^ t K . . .,w in ) 

%<] J 

where Py ln is the Macdonald polynomial with parameters tu = q~ 2 and qu = t~tjf ■ 
The proportionality factor is the leading coefficient, i.e. 

li,n = /^M*) ([!]!)" 

Proof. We have proven above that the correlation function in the theorem is a homogeneous 
symmetric polynomial that satisfies the wheel condition, more precisely the (2, £)-wheel con- 
dition. Therefore it lives in a vector subspace spanned by Macdonald polynomials associated 
with (2, £)- admissible Young diagrams. 

The polynomial depends in 2n variables and has a combined degree of inin — 1). The 
diagram Y^ n is the smallest one that satisfies the wheel condition (with 2n parts) and it has 
exacly £n(n — 1). Therefore, it corresponds to the only Macdonald polynomial with such 
conditions. □ 



As (k\e PJ-\k) also satisfies the wheel condition, we can express it as well as a sum 
over Macdonald polynomials: Wi) (k\ e 2to/3 J r \k) = c^yPy^i, . . . , wg n ) where the 
sum runs over all (2, £)-admissible Young diagrams with in(n — 1) + n(£ — k) boxes. By 

y(*0 

proposition 3.1, we know that the leading monomial is w e < n and therefore only term appears 
in this sum: 

' In \ 

JJwA (k\e 2en ^\k) = 1 ^P y(k) (w 1 ,...,w en ) 



v i=l 



with coefficient 



Notice that this implies that P y (o) = (JT, w i) Py e 



3.1.2. Derivation of integral formulae. We have computed all the ingredients to compute 
fen' ^ n orc ^ er ^° express the result, we use the following reparameterization of the 
index sequence {bi, . . . , &2n}- Let e = {ei, . . . , ee n } be a non-decreasing sequence, such that 
1 < Q < 2n, and such that, the same value can not be repeated more that I times. Then, 
there is a bijection between the set of neutral spin sequences {bi, . . . ,b 2n } and the non- 
decreasing sequences e, given by: bj = #{ej such that Cj = j}. To each e« equal to j 
corresponds a current f{wj) that is used to raise the spin of <&o(zj). 

We then compute: 
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• From the product $0(^1) • • • &o( z 2n) we get: 



%( Zl ) . . . $ (> 2 „) = TT(-Zi) 2 ~A ! Z_2° • $ ( 2l ) . . . <!> ( Z2n ) : 

J- I n O \ n ZA a \ 



i<j \q Zt , q 

• When ei < j, we have pairs of the form f(wi)§o(zj), and then: 

/(O$ (%) = ^-^:/K)$ (^): 

• When €i > j, we have pairs of the form §°(zj)f(wi), and then: 

Moffat) = -jrrz ; : /K) $ o(gj): 

• Finally, we get an extra term for the case 6j = j : 

ttiW _,/.../(n^^)*ofe)n/( 

where the last product is in the increasing order. 
We conclude using lemma 2.1. We first write the case = 0: 



W; 



r (°) TT. TT ^§ W ■ 



(3-2) <L 62n (^i,---,^)=7^n^ n 



t=l l<i<i<2n |^ +2 2-;? 4 

T^r WjrfWj T-r U>j - U>j Py itn (Wi, Wj n ) 

f = [ 27^^ f<j W l ~ ^ IX^i - rii> ej K ~ 1 lz i) 

Where the countours are chosen such that 

\u)i\ < \q~ £ Zj\ for all j < 6j |iUj| > \q Zj\ for all j > 

> |g 2 -2j| for all z < j \wi\ > \q~ 2 Wj\ for all i < j 

This choice is consistent with remark 2.4. 

The same computation can be repeated for = (&| . . . $b 2n and we obtain: 



2n 



(3.3) . . . , = tS n^(-^) f n 



9 



00 



X 



1 / / 

11 97T?' 11 7/1, — 



i 2tu 1 Wj . _ q 2 w . ( q i Wi - Zj ) Yl.^ ( Wi - q i Zj ) 



In general, we do not expect a simple formula for the Macdonald polynomial -P v o) , except 
in two cases: 
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(1) If i = 1, the parafermionic part of the current is trivial, and P Y (k) = Yli w l~ k 

Yli<j( w i ~~ <l 2w j)( w i — Q~ 2 Wj), and we recover the standard level 1 formulae found 
in e.g. [17]. 

(2) If i = 2, the parafermionic field reduces to a (g-deformed) fermionic field [3], and 

P <; = T^ggf* Pt ( (S^gg^) ),^. "tee is 
the elementary symmetric polynomial, reproducing results of [15]. 

In the rest of this section, we shall write \I/ := when there is no risk of confusion. 

3.2. Quantum Knizhnik— Zamolodchikov equation. A detailed proof of the qKZ equa- 
tion can be found in [16, appendix A]). For the convenience of the reader, we provide a short 
graphical proof in appendix B.2. 

Define 

R±(z) = 4>±{z)R{z) 

with 

6 ( Z ) = ^ 2 / 2 (g 2 *; ^ 

0+(z) = MV*)" 1 

Note that compared to the other normalization R(z) we have used before, one has: 

f (-l) m £ = 2m even 

[(-lHz/g) 1 / 2 v (,^4)^4/,.^)^°° ^ = 2m+lodd 

and 4>+(z) = (p-il/z)- 1 . 

The following equation is then satisfied by ^{z\, . . . , z£) = (k\ $(^i) • • • ^(^x) \k): 

(3.4) i?_(z i+1 /(s^)) . . . R_(z L /(sz i ))Kl +k R+(z 1 /z i ) . . . R + (^-i/ ...,z L ) 

= fy(zi,...,szi,...,z L ) i = l,...,L 

where s = g~ 2 ( £ + 2 ); all the indices are suppressed and can be recovered by following the spec- 
tral parameters Zi, see appendix B.2 for details. This is the quantum Knizhnik-Zamolodchikov 
equation, first introduced in [13]. 

As a consequence of (2.15) for the VOs we also have the following identity (exchange 
relation): 

(3.5) #M+i(zi/*i+i)*(si, ...,zl) = • • • , z i+i, Zi,..., z L ) 
We shall come back to this equation in section 3.4. 

3.3. Recurrence relations. Define for convenience the "dual vertex operator" $*, which 
by self-duality of evaluation representation can be expressed in terms of $ as: 

(3.6) $* b (z) = Ci - b $ t - b (q- 2 z) 
where c b = (-l) 6 ^ 6 "'" 1 ) [&]![€ -b}\/[£}\. 
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From general arguments one can prove that 
(3.7a) = Vtf* 

(3.7b) = A* 

6=0 

where is a scalar which depends only on the highest weight representation, and whose 
explicit expression, given in [ , Prop. 4.1], will not be needed in what follows. 

This implies immediately the following recurrence relations for the correlation functions: 
(3.8a) M'...,, v ....(. ..,z, q~ 2 z } . . .) = ^S bW/ *...(. . .) 

£ 

(3.8b) d-b^...,b,e-b,.X- ■ ■ , 9" 2 ^, z,...)=g k *...(. . .) 

6=0 

where . . . are omitted (arbitrary) arguments. 

3.4. Polynomiality, cyclicity. It is convenient to redefine 

(3.9) ^\z 1 ,...,z L ) = f[4 L - im+k/2 II F(z j /z i )^ k \z 1) ... 1 z L ) 

i=l l<i<j<L 

where 

= ^/ 4 



(g 4 z; g 4 ) 



Note that the prefactor Ylt=i z^ L ~ 2 combined with equation (2.14) implies that 

*l?(xzi, • • • > %Zl) = x in ( n ~ 1 ) +kn Sfr (zi, . . . , z£) (where L = 2n); in fact, rewriting the integral 
formula (3.2) as: 

2n £ „ „ £n j 

*)-£n<M? n n<«***-*>f-/nig 

i=l \<i<j<2nr=l J J i=l 

P (k){W 1 , . . .,W en ) 



n 



Wj — Wi 



(where we recall that e = (ei, . . . , e&J is the sequence such that bj = #{ej such that £j = j}) 
one can show (see appendix D) that is a (vector-valued) polynomial in the variables 
Zi, . . . , Zl of degree £n(n — 1) + kn. Furthermore, its coefficients are rational functions of q 
whose denominators are products of 1 — g 2j , j — 1, . . . , I. 

For example, when I = 2 and n = 2 there are nineteen components, among which: 
*i°2,i,o0i> Z2, z s , Za) = {q + q 3 ){zi - q 2 z 2 )(z 2 - q 2 z 3 )(z 3 - q 2 z 4 )(z 1 - q 6 z 4 ) 
*i°2,o,i(>i> z 2, z 3 , za) = -(q + q 3 )(zi - q 2 z 2 )(z 3 - q 2 z A ) 

x {q 2 z x z 2 + z x z 3 - q 4 Z!Z 3 - q 4 z 2 z 3 - q 4 z 1 z 4 + q W z 3 z A ) 
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For generic n and I the simplest component is the one with hi = for all i < n and hi = I 
for the remainings % > £: 

n £ 

(3.H) *t,oA...s = (-iy nj ^ q ^ n+i ^ n * n ik^ - - ^ 

i=l 0<i<j<n r=l 

and all its rotations. 

Now define yet another normalization of the R- matrix, namely R(z) = -J^j^R(z), which 
is adapted to St from (3.9). Explicitly, 

The exchange relation (3.5) becomes 

(3.13) Ri,i+i{zi/z i+ i)^(zi, . . . , z L ) = &(zi,...,z i+ i,Zi,...,z L ) 

Starting from (3.4) at say i — 1, switching to St, applying repeatedly (3.13) and using 
4> + (l/ z)F(sz) / F(l/ z) = (— 1) , we find the cyclicity relation: 

(3.14) a*,..., ZL _ X ) = 9(zi, z 2 ,..., sz L ) 

where S : V ZL <g> V Zx <g> • • • <g) y^ L1 — >• <8> • • • <8> V SZL rotates cyclically the tensor product. 

The system of equations (3.13-3.14) is similar to the one that appeared first in [28] in 
the study of form factors. Note that the power of s is a reflection of the homogeneity of 
*&(zi, . . . , Zl). In components, (3.14) writes: 

( _ 1) ^(„-lW2 g (l +fe )(2 6i ^)^ )6i] >bL _ i{zu Zu . . . , ZL _ l} = * bl _ bL ( Zl , Z 2 ,..., SZ L ) 

Other formulae can be rewritten in terms of \I/ as well; in particular, the recurrence 
relations (3.8) become: 

(3.15a) *.„,6,y > ...(- Zi+i = q~ 2 Zi, ■■■) = 

i-l I LI 

9k n n& - n - <r 2(r+i) ^)* .(• ■ o 

j=Xr=X j=i+2 r=l 

t 

(3.15b) ^ c e-b^...,e-b,b,-(- ■ ■ > z % = q~ 2 z i+1 , z i+1 , . . .) = 

6=0 

i— 1 t LI 

k n - n - ?- 2(r+i w)* .c . o 

j=l r=l j=i+2r=l 
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3.5. Wheel condition. Finally, we point out that satisfies a wheel condition of a different 
nature than that of the Macdonald polynomials of sect. 3.1.1. Namely, we have the following 
theorem: 

Theorem 3.4. Assume that three parameters z^, Zi 2! Zi 3 , i\ < %2 < 13, form a "wheel": 
z% = q 2<fiJrb >z, Zj = q 2b z, Zk = z with a, b integers such tht a, b > and a + b < £ + 2. Then 

(3.16) *(..., q 2(a+b) z, q 2b z, ...,z,...) = 

The proof is similar to the one given in [32, theorem 3] and is reproduced in appendix C.3. 

Note furthermore that the entries of ^ are nonsymmetric polynomials, hence the ordering 
condition for its arguments in the wheel condition. This wheel condition is not a special case 
of the one considered in [18] for nonsymmetric Macdonald polynomials. 

4. Application to integrable spin chains 

In all the discussion that precedes concerning vertex operators, it has assumed that \q\ < 1. 
However, we have seen in section 3.4 that one can get rid of all infinite products in explicit 
formulae by a redefinition of the solution of qKZ and of the i?-matrix R. The former 
becomes a polynomial of Zi,...,Zl (with coefficients of the form P(q)/Yii=i(^ ~ q 2t ) ki ) 
whereas the latter becomes a rational function of them. As a result, we can now relax the 
constraint \q\ < 1 and in particular consider the case where q 2 is a primitive (£ + 2) th root 
of unity, so that s = q~ 2( - e+2 ^ = 1. It is the purpose of this section to show that SI/ then 
becomes an eigenvector of an integrable transfer matrix. 

Note that since SP ~ (k\ $ . . . $ \k), our procedure is similar to a "matrix product Ansatz", 
where the role of the matrix algebra is played here by the Zamolodchikov-Faddeev algebra 
50, 10, 21]. In fact, a similar procedure has already been proposed in [1, 2, 19]. However, 
there is a crucial difference between our situation and theirs. In the aforementioned papers, 
in order to impose periodic boundary conditions, a trace is taken, which would correspond to 
computing tr($ ...$). Unfortunately such a trace is divergent in our setting. Also observe 
that such an Ansatz would produce eigenvectors of integrable models for an arbitrary value 
of the quantum parameter q, a claim which we do not make here. Instead we take a vacuum 
expectation value of VOs; the price to pay is that rotational invariance (periodic boundary 
conditions) is only restored at special roots of unity, namely, q 2( ~ t+2 ^ = 1. One can also 
regularize the trace by adding a x d : tr(x d $ •••$), producing finite temperature correlation 
functions, as already mentioned in the introduction; but that would also spoil the rotational 
invariance, and is therefore not directly relevant in the present context. 

4.1. The model. Define the inhomogeneous monodromy matrix to be an operator on V z £g> 
V Z1 <g) • • • ® V ZL of the form: 

T(z) = RioCzi/z) • • .R L o(zl/z) 

where the index corresponds to V z and non zero indices to V Zi , and the i?-matrix is given 
as before by (3.12). with R = "PR. 
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The twisted transfer matrix is then defined by taking the trace over the auxiliary space 
V z , with a twist (in which it is convenient to absorb the sign (— 1)^): 

(4.1) T(^)=tr ((-g 1+fc rT(^)) 

The Yang-Baxter equation implies that transfer matrices commute for different spectral 
parameters: 

[T(,),T(z')]=0 

The diagonalization of T(z) is usually performed using Bethe Ansatz, which produces 
eigenvectors defined in terms of "Bethe roots", i.e., solutions of certain algebraic equations. 
We do not pursue this route here. 

A local Hamiltonian can be extracted from the transfer matrix in the homogeneous limit 
where all the Z{ coincide. Assume that Zi = 1. According to (3.12), R(l) = V, and therefore 

(4.2) T(l) = (-q 1+k ) a °' L S 

where we recall that S is cyclic permutation of factors of the tensor product. T(l) is a discrete 
analogue of the momentum operator. Expanding to next order, we obtain the Hamiltonian: 



where h jd+1 = \£R j>j+1 (z)\ z=1 , i < L, and h L>1 = \(-q l+k Y aa < L ± &m(z)\ z=1 (-q l+k y^ 
(twisted periodic boundary conditions). The factor of 1/i is introduced for convenience, see 
below. 

4.2. The simple eigenvalue. Assume that q 2 is primitive (£ + 2) th root of unity (primitive- 
ness being necessary for ^ to be well-defined). Consider given by (3.10). An important 
remark is that the values of qu = q 2 ^ +1 ^ and tu = q~ 2 coincide, so that the Macdonald 
polynomial appearing in this expression becomes simply the corresponding Schur polyno- 
mial. Therefore, 

2n I „ „ in , 

•£*.(*..... w-«n«Mf n n(«**-*>/-/n^ 

^4 3^ i=l l<i<j<2nr=l J J i=l 



X 



Z 3, 



where hj = Lfj + L^J + 3- 

Next, note that (3.14) can be simplified if one assumes s = g~ 2 ^+ 2 ) = 1 to 

(4.4) (-q)( 1+k ^S*(z 2 , ...,z L , z x ) = ...,z L ) 

Now consider the effect of the transfer matrix on the correlation function *f?(zi, . . . , zl)- 
Writing the transfer matrix T(z) as the trace of the monodromy matrix T(z) and using 
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(3.15a), we have 3 

(4.5) T{z)*{ Zl , ...,z L ) = IT 1 ]T((-g 1+fc rVT fe ,^) — *m-6V..(^, q~\ z x ,...,z L ) 

where II = firfe Il}=i Yl r =i( z j ~ <? 2r ^) Hj=i+2Hr=i( z j ~ g~ 2(r+1) ^) is the factor appearing in 
(3.15a), and the ... in subscript mean that only the first and the last index are fixed, the 
rest forming a vector in V Zl <g) • • • <S> V ZL . 

Using (4.4), we can rewrite this as 

T(z)^f( Zl , . . . , z L ) = ir 1 V T b , b/ (z)—* £ _ b ,_ b (q- 2 z, z u ...,z L ,z) 

TP Ce ~ b 

Finally, writing T{z) as a product of R- matrices and applying the exchange relation (3.13) 
repeatedly, we find: 

T(z)¥(zi, . . . , z L ) = IT 1 V —tyt_ b>K ..(q- 2 z, z,zi,...,z L ) 

h 

= *(z u ...,z L ) 

where the last equality follows from (3.15b). 

Noting that ^ ^ because of (3.11), we conclude that if q 2 is a primitive (£ + 2) th root of 
unity, *&(zi, . . . , z£) is an eigenvector of the (inhomogeneous, twisted) transfer matrix T(z), 
with a trivial eigenvalue. A graphical interpretation of this proof can be found in appendix 
B.3. 

The twist is — q 1+k , but note that only its square is meaningful. As k varies from to £, 
q2(i+k) S p ans all £ + 1 nontrivial (£ + 2) th roots of unity. 

We conjecture that if q = — e ±l7r /^+ 2 ) (for all twists — q 1+k ), the eigenvector we have just 
constructed corresponds to the largest eigenvalue of T(z) for z,- L of modulus 1 and sufficiently 
close to 1. In particular, for Zi = 1, we conjecture that ^ is the ground state eigenvector of 
the Hamiltonian H (with zero ground state energy). For a discussion of the latter statement, 
see sect. 4.4. Note that at q — ± e ±l7T /( e + 2 ) ; the spectrum of H is real. 

4.3. Relation to loop model. In [ ], a "higher spin" loop model was built by fusing the 
standard Temperley-Lieb loop model. Since the latter is related to the spin 1/2 represen- 
tation of U — U g (sl(2)) in the sense that it is equivalent to the XXZ spin chain/6- vertex 
model, the fused loop model must be related to higher spin integrable Hamiltonians/transfer 
matrices. The only nontrivial issue is that of boundary conditions. Let us discuss this here, 
first in the spin 1/2 case, then in the fused case. 

4.3.1. Equivalence to spin model and twisted boundary conditions. We first define a link 
pattern of size L to be a planar pairing of {1, . . . , L} viewed as boundary points of a disk, 



'The argument that follows was suggested to us by R. Weston, to whom we are indebted. 
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e.g., 



32 32 32 32 32 




> 



56 56 56 56 56 

V / 



We shall identify a link pattern with the fixed-point-free involution which sends paired points 
to each other. Note that L has to be even for Cl to be nonempty. 

The span of link patterns, C[£l], can be identified with a subspace of (C 2 )' X)L ; indeed, 
there is a linear map ip which to each 71 G Cl associates a vector in (C 2 ) <Xli , which is given 
by the following explicit formula: 

(4.6) E II <w (I) ,i(-<?r i+ k 

€€{0,1} L l<i<vr(i)<L 

where on the r.h.s. the v e are the standard basis of (C 2 )® L . In appendix C, we show that tp is 
injective, and that for generic q it is in fact an isomorphism from C[£l] to the L^-invariant 
subspace of (C 2 ) (g)L (recall that U\ is the horizontal subalgebra of U). 

The change of basis described by (4.6) has the following diagrammatic interpretation: 
one sums over all possible orientations of the arcs, assigns a weight which is equal to 
2^ total angle spanned by the arcs ^ an( ^ reconstructs the spin state by recording the ori- 

entations at the endpoints of the arcs, identifying v\ = v-f, vq = v±. 

The actual shape of the disk and the positioning of the points on its boundary are irrelevant 
and result in gauge transformations (i.e., moving around the twist). The choice above 
corresponds to all points aligned on a straight region of the boundary, e.g., 



2 




1234 1234 1234 1234 



— — q Vf^u + v nn + v im ~ Q v un 

The simplest way to compute the twist is to consider the rotation one-step to the left of 
link patterns, S, and map it to the spin space. We find that S = (—q) a °- L S where S is 
(untwisted) spin rotation, compare with (4.2) at k — 0. 

For the sake of completeness we also indicate how to map covectors (and therefore, opera- 
tors) . It is convenient to think of covectors in the loop picture as planar pairings outside the 
disk. 1 The pairing between vectors and covectors is to paste them together, and assign them 

4 Note the asymmetry: the region outside the disk is homeomorphic to an interval times a circle ( "identified 
connectivities"), not a disk ("unidentified connectivities"). It would be tempting to add the point at infinity 
to restore the symmetry, but that is not allowed because we want the winding angle of any closed curve to 
be 2?r. 



20 



T. FONSECA AND P. ZINN- JUSTIN 



r to the power the number of closed loops thus formed, where r = — (q + q 1 ). Note that r is 
nothing but^^+^^, where we have used the same rule as before that the weight of a line 
is equal to (— g)^F total angle s P anned b y the arcs . Therefore, the exact same graphical rule must be 
used for covectors: sum over all orientations, give a weight of (— q) 2^ total anglc s P anned b y the arcs 
and build the dual spin state from the orientations at the endpoints of the arcs, with 



Sij. 



4.3.2. Fusion. We shall not reproduce here the fusion procedure, i.e., consider representation 
theory of U, but rather consider only the non-affine part U\. We therefore view the C/i-module 
V = <C e+1 as a submodule of (C 2 )®^, as well as the projection p : (C 2 )®^ — > V (see C.2 for 
its explicit definition). 

The £-fused link pattern of size L of [ ] is then obtained from ordinary link patterns of 
size £L by grouping together vertices into groups indexed by = 1+ L^y-J an d by restricting 
to link patterns with no connections within a group, i.e., 



£t,L = {vr G L lL : 7r(i) = j ^r,i^ r,} 



For example, 



C 



3,4 




1 • 




In the fused case, to a 7r E C^l is then associated the corresponding vector after projection: 

L 

(4.7) ^(tt) = Y[pi(p{ir) 



where pi is p acting on the factors £i + 1, . . . , i{% + 1) of the tensor product (C 2 )® iL . 

Clearly, twist and fusion commute, so that we have the same twist (— 1) K\ = (—q) a ° in 
the fused case. 

In appendix C, we show that the map tpg is well-defined and injective provided q 2 is not a 
r th root of unity, 2 < r < £, and that for generic q it is in fact an isomorphism from C[£^ ; l] 
to the t^i-invariant subspace of (C 2 )® L . 

The integrable inhomogeneous transfer matrix can be formulated directly in the language 
of loops, as is done in [32]. Once converted to the spin space, it becomes exactly the twisted 
transfer matrix T(z) discussed in 4.1, with twist given by k = 0. 

4.3.3. Some results. Consider now ^°\zx, . . . , z£). It is proportional to \I/^ ^ (^i , . . . , z£) = 
(0| • • • $(zl) |0), and therefore, from the intertwining property of &(z), is U\ -invariant. 

We conclude that it is in the image of ipi for generic q, and by continuity it is so for q 2 
(£ + 2) th root of unity. We can therefore consider its preimage tp t By construction 

it is an eigenvector of the transfer matrix of the loop model, with eigenvalue 1. And it is a 
polynomial of degree £n(n — 1) in the Z{. At the particular value q = —e~ ln ^ i+2 \ which is 
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the only one considered in [ ], coincides with the eigenvector studied in that paper, 
and its degree is the content of Conjecture 1 of [32], which is thus proved. 

4.4. Relation to supersymmetric fermions with exclusion. In [ ], a supersymmetric 
model of fermions on the one- dimensional lattice was introduced. It has the exclusion rule 
that at most £ consecutive sites may be occupied. In the case £ = 1, it was found to be 
equivalent to the XXZ spin chain at A = —1/2 [29]. More generally, in a recent paper [14], 
Hagendorf conjectured that for any £, this fermionic model was equivalent to the integrable 
spin 1/2 chain at q = e ln ^ e+2 \ Note that models related by q — > —q and q — > 1/q are 
closely connected (in particular the choice q — > —q depends on the sign convention for the 
Hamiltonian) . It was also stated that the admissible twists are of the form e 27rim /i i + 2 ) a o ; 
where m = 0, 1, . . . ,p, p = £+ 1 for £ odd and p = £/2 for £ even. However, it seems difficult 
to compare his twists with ours because the equivalence is only conjectural and not explicit 
in the general case, and may modify the twist. See also the computation of the Witten index 
in [14, sect. 3.5]. 

Adapting this discussion to our setting, we conclude (conjecturally) that the Hamiltonian 
H at q = —e ±lw ^ i+2 ^ should be supersymmetric for all twists of the form ±q n . Furthermore, 
it should have a supersymmetric ground state exactly when n^O (mod £ + 2). This would 
imply that our state which has zero energy, is indeed the (supersymmetric) ground state 



Here we follow the procedure of Ding and Feigin [8], we intend to construct a level £ Drin- 

feld realization of U q (sl(2)), and to rewrite the current in terms of bosonic and parafermionic 
operators (these being appropriate g-deformations of Conformal Field Theory operators, cf 
the parafermionic fields of [31]). 



of H. 



Appendix A. Construction of (/-deformed parafermions 



Let 



E^ 1 /(*) 



z 



(Al) 



k + (z) 



E k t^ k ~( z ) 



—i 



z 



We want to construct operators that satisfy the following relations: 



g{q c z/w) 

(A.2) k ± (z)e(w) = g (q- c/2 (w/zf 1 )^ e(w)k ± (z) 




k ± (z)f(w) = g {ql* (w/zf'f 1 f 



e(z)e(w) (z — q 2 w) = e(w)e(z) [q 2 z — w) 
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f( z )f( w ) {(fz ~w) = f(w)f(z) (z - q 2 w) 
(q - q-^zwfciz), f{w)} = 5 (q~ c z/w) k + (q- c/2 z) - 5 {q c z/w) k~{q c/2 z) 
where c is a central element. The functions g(z) and S(z) are defined by: 

Remark A.l. The Chevalley generators are given by: 

K = q c K E = f l (k+)- 1 F = k+e-i 

Ki = k+ E 1 = e Fi — f 

A Hopf algebra structure is obtained as follows. We define the co-multiplication: 

A(q h ) = q h ®q h 
A(e(z)) = e(z) <g> 1 + k-(q Cl/2 z) <g> e(q Cl z) 
(A ' 3) £(/(*)) = f(q c *z) ® k + (q^z) + 1 ® /(z) 

A(A; ± (^)) = k^q^z) ® k^q^'z) 
where C\ — c <S> 1 and C2 = 1 <S> c; the co-unit: 

(A.4) e(q h ) = 1 = i(k±(z)) e(e(z)) = = e(f(z)) 

and the ant ip ode: 

~a{q h ) = q- h 
a(e(z)) = -k~{(f ct2 zy x e{(f c z) 
~a{f{z)) = -f(q~c)k+(q- c / 2 z)- 1 



(A.5) 



A.l. Highest weight modules. We define a highest weight module as the unique module 

V(A) := U q (sl(2))v\, where A = n A + niA x is the highest weight, and v A is the highest 
weight vector, which satisfies: 

EiV K = K iVA = q^v A F t (A ' a ' )+1 v A = 

the inner product being defined by (Aj,<x,-) = 5 it j. Except in the trivial case A = 0, the 
modules thus generated are irreducible and infinite-dimensional. 

We will give a more explicit way of constructing these modules depending on the eigenvalue 
£ = n + n\ of the central element c. 

A. 2. Level 1 realization. We first build a realization of the U q (sl(2)) algebra for the case 
c = £ = 1 in terms of bosonic currents. We introduce a new set of bosonic operators, which 
satisfy the relations: 

r 1 _ X \M M F ] nr 

[a n ,a m \ — o„+ m ,o [a n ,0!\ — Ao n $ 
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Then we can build the currents by: 



, , „-"/2 a -n -n V/,-°/ 2a » 



, , _"/2 n -n _n v»/2%,-n 

k + (z) = q ao e {q - q ~ 1) ^ a " z ~ n 

k~(z) = g" a o e - (<? - 9 " 1)Ea -" 2 " 

these satisfy all the above relations. We omit such computations. 

They act on the usual bosonic Fock space %x = © AgZ ((llm>o a -m ^y the following 
rules: 

<hf \h) B = aj \h) B if i < e a f \h) B = f \h + 2) B 

&if \ h ) B = \ a ii f\\ h ) b if*>0 oof \h) B = fh\h) B 

where / G C[a_i, a_ 2 , . . .]. 

This produce two highest weight modules, depending on the parity of h, corresponding to 
the two highest weight vectors 

^A = |0) B V Al = \1) B 

A. 3. Level £ realization. Let us denote ax : U — > C(Hx) the level 1 realization. We 
produce, using the coproduct A, a level £ representation as the tensor product of £ copies of 
V, i.e., if we define A e = (1 ® . . . <g> A) o . . . o (l (g) A) o A, then a = A^(o"i CS> • ■ • ® ax) is the 
representation on H e = Hf e . 

We also define currents in Tie in the obvious way: 

e(*) = £<(e(*)) f{z) = K\f{z)) 
k±{z)=h i {k ± {z)) q c = K\q c )=q c ®...®q c 

As A is a co- multiplication, this will satisfy automatically the relations (A. 2), with c = £. 
In that way we can compute e(z) for any £. For example, if £ = 2 we get: 

(A.8) e(z) = e(z) g> 1 + k-(q Cl/2 z) ® e(g C1 ^) = e 1 ^) + e 2 (^). 

where is the value of c in V^. 

We perform a mode expansion of e(z), for generic £: 

(A.9) e(z)=5>*(z) 
where (using q = 1) 

(A.10) e\z) = (g) k-tf-^z) <g> e^z) (g) 1 

where j indicates the position on the tensor product. 
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We repeat the process for f(z): 
(A.ll) /(*)=£/*(*) 

where (using q = 1) 

(A.12) f(z) = (g) 1 ® /(^) (g) fc+(^+^) 

The same goes for fe^z): 

(A.13) fc ± (z) = (g)fc ± ( 9 ±1 ^ ±i z) 

Let us denote |fci, . . . , fa) — \ki) B <g>. . ■ <S>\kg) B , where fa = 0, 1. Let k = ^ |&i, • • • , fa), 
as a tensor product of highest weight vectors, is a highest weight vector with highest weight 
A = (£ — k)A + kAx. Note that for a given A, except for k = or I, this construction 
is not unique. For our purposes, any such highest weight vector provides a realization of 
the highest weight module, and we will use the notation \k) for any of them (our results 
do not depend on the choice, i.e., when we permute fa -B- kj). Moreover, these vectors are 
normalized such that (k\k) = 1. 

A. 4. Parafermions. The computation of commutation relations between e l (z) (or f l {z)) 
and fe ± (w) it is not complicated, although it is long. The result for e l {z) is: 

e i (z)k + (w) =: e i (z)k + (w) : 



(A.14) 



and for f l (z) is: 



(A.15) 



e\z)k-{ W ) = z ;*JZl ■ ei ww ■ 

k-{w)e i {z) = q~ 2 : e\z)k~{w) : 



f(z)k + (w) =: f(z)k + (w) : 
k + (w)f(z) = ^ q _ : f\z)k+{ W ) : 

z _ n 2 + e / 2 W 

r(z)k-(w) = z _ q q _ 2+e/2w :f{z)k-{ W ): 
k-(w)f i (z)=q- 2 :f(z)k-(w): 



Surprisingly, these results do not depend on i. This fact, they tell us that we can factor 
e l (z) and f l (z) into two parts: one (parafermionic part) that commutes with ^(w) and 
depends on i and a second part (bosonic part) that does not depend on i and does not 
commute with fc ± (w). 
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Following Ding and Feigin, let us build the bosonic part. Define (let n be a non-negative 
integer) : 



a±n = a ±n 



< n = 01® q~ n ^a ±n (g) 1 

j<i j>i 

(correcting a sign misprint in [8]). 

We could think that a n is only the result of applying the co- multiplication to a n . But this 
is false. In fact, there is no Hopf algebra for these operators. 

These new operators satisfy the following commutation realation: 

(A- 1 ' J [ a ri; a mJ — Cn+m,0 

It can be easily proved. 

The main reason for this definition is the fact that: 
(A.18) fe ± (z) = g ±a o e ±(g-< ? - 1 )E„-±^" 

where, for example, q a ° means q a ° ® • • • <8> q a " ■ Notice that this is the natural generalization 
ofife^z). 

Now the idea is to factor e t (z) and f l (z) using these new bosonic operators. That is, we 
want to write them as e l (z) = t(z)^ +t (z) and f l {z) = f(z)£ _i (z), where t(z) and f(z) are 
expressed in terms of a± n and e a = . e a (compared to the notations of the main text, we 
have renormalized the zero mode as ex. = 2£{3). 

The following Ansatz is made: 

(A.19) it 

In this way, the commutation relations of t(z) (or f(z)) and fc ± (w) have the form of (A. 14) 
and (A. 15), i.e., 



0= [£+*(*), fc±(w)] 
(A.20) . ^ J 

0=[rW,fc ± H] 

for all i. 

In fact, one can prove that 

(A.21) [a n ,£ H (w)) = [a,£ H (w)] = 

for any integer n and i G {1, 2, . . . , £}. 

The operators £ + (z) = X^C 4 " 1 ^) an d £,~( z ) — J2i£,~ l ( z ) are the so-called parafermionic 
operators. 
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We can split a vector \k) into a bosonic component and a parafermionic component. The 
first being an average \k) B = \ k /f) B <8> • • • <8> \ k /f) B - This requires that we enlarge the Fock 
space to Hi = hgz/ , <(n m >o a -m I^)b>- 

A. 5. The wheel condition. These operators satisfy an important relation. 
Theorem A. 2. e(z),f(z) satisfy the wheel condition: 

e(z)e(q 2 z)...e(q 2£ z) = 
f{q 2 "z) . . . f{q 2 z)f{z) = 

This is theorem 2.5 in [8]. 

In order to prove this fact, we need the folowing relations for e l (z) (let j > i): 



(A.22) 



_2 

e\z)eHw) = z ~ q w : e i (z)e j (w) 
z — q 2 w 



e i (z)e i (w) = q 2i 2 (z - w)(z - q 2 w) : e'^e^w) : 
e j (z)e i (w) = q~ 2 : e j (z)e i (w) : 
And the equivalent for f l {z): 

f(z)f(w) = Z -f^ : f(z)f(w) : 
q A z — w 

(A ' 23) f(z)f(w) = q 2l ~ 2 \z -w)(z- q 2 w) : f(z)f(w) : 

f(z)f(w) =: f{z)f{w) : 

Note that when u> = g 2 z (or u> = q~ 2 z for the /*(-z) case) these expressions simplify. This 
is our main tool to prove theorem A. 2: 

Proof. The product e(z\) . . . e(z n ) can be written as the sum over all possible decompositions: 

e( Zl ) . . . e(z n ) = J2 e ei ( Zl )...e en (z n ) 

Ei,...,e„ 

The regular part, that is the : . . . : part, has no poles (aside and oo). Thus, setting 
w = q 2 z, we see that only e^{z)e l {q 2 z) survives (we are assuming j > i). 



2n-2 z) 



e(z)...e(q 2n - 2 z)= e ei (z)...e e 

ei>...>e„ 

but such sequence t\ > . . . > e n is impossible if n — £ + 1, considering that c t G {1, . . . , £}. 
The result follows. The proof for f l (z) is identical. □ 



Appendix B. Graphical calculus 

We describe here the graphical calculus to represent U q (sl(2)) invariants. It is a convenient 
way to derive various relations while keeping track of the spaces involved, i.e., it acts in a 
similar way as a "type checking" tool for programming. 
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B.l. Basics. Vector spaces are represented by oriented lines, such that tensor products 
correspond to stacking lines from left to right if looking in the direction of the orientation. 
Highest weight representations are denoted by thick lines, evaluation representations V z 
by thin lines with a label z. Thus, a (type I) vertex operator &(z) becomes: 



and its (vacuum to vacuum) correlation function is: 



= (0|$(zi)...$(z L ) |0) =o 



Z\ . 


zl . 




) 




-4— ( 



where the white dots on the left (resp. right) correspond to (0| (resp. |0)). 

We need two more objects obtained by applying representations to the universal i?-matrix 
K (or to V(U- 1 )): the i?-matrix 



R+(z 1 /z 2 ) 



/z\ 



R-{z 1 /z 2 ) 



Z2 S 



where the argument is the ratio of spectral parameters of left and right incoming lines (note 
that R + and i?_ only differ by a normalization); and the L- matrix 



-1/2, 



L+(z) 



L_(z) = 



. q^z 



q e / 2 z 



-t/2. 



a(L + (z)) 



a(L_(z)) = 



q*/ 2 z 



-1/2, 



-111, 



1 

q e / 2 z 



Note that thin lines pick up a factor of q ±E when they cross thick lines. 

Lines can also be slid across intersections of other lines (Reidemeister move III) due to 
the Yang-Baxter equation, or across the trivalent vertex of a VO (L/VO commutation). 

Finally, it is easy to check on generators that the square of the antipode satisfies a 2 (x) = 
q-ao+4d X qao-'id f or a jj x j n jj 'j^jg j m plies that any line that does a full ±27r rotation can 
be replaced with a straight line with an insertion of q ± ( a o- id ) times a central element, e.g., 



(B.l) 



P 



c <> 



C -+ 



where the dot represents insertion of g ao_4d . We shall not need the values c, C of the cen- 
tral element in what follows. The q~ 4d acts as a multiplicative shift of q~ A of the spectral 
parameter, which propagates along the lines. 
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B.2. Proof of gKZ equation. We begin with a local relation: consider the following equal- 
ity, obvious diagrammatically: 



\q- 2t z 



q -z 





■ —IV 
q 












If we use (B.l), recognizing conjugation by q a ° Ad on the thick line, the coproduct A(g a ° Ad ) 
qa -4d gao-id an( j g na Uy fa e intertwining property of the VO, we can simplify this to: 







sz 


. J 







. sz 




— (- 



where s = q 2 ( i+2 \ and we recall that the dot represents insertion of K\ = q a ° accompanied 
by a shift of the spectral parameter z — > q~ 4 z, which we pulled out of the thin line. 

Now apply this relation to the correlation function {z\, . . . , z£) described graphically 
above, and use the highest /lowest weight property of |0) and (0|, i.e., L + (z) |0) = 1 and 
<0|L_(*) = 1: 



(B.2) 



o 



zi , 


SZi . 


Z L . 




)— 




<- 


c 



o 



1 



L 



-o 



which is exactly the quantum Knizhnik-Zamolodchikov equation (3.4). 

In case (0| . . . |0) is replaced with (k\ . . . \k), k = 0, . . . , £, there is an additional contribution 
coming from the zero modes of L±(z) ~ g a o®a-o/2^ go ^ e - m fa e ,-[ ^ becomes g( 1+fc ) a °. 

The qKZ equation is valid for the correlation function of a product of arbitrary VOs; 
however, in the special case of perfect VOs, we have the additional relation (2.15): 












— «- 





where the flat crossing corresponds to another normalization of the i?-matrix, denoted by 
R(z) in the text. This in turn implies the exchange relation (3.5) for correlation functions, 
and, combined with the qKZ equation at i = L, produces the cyclicity property: 



o 





zl-i , 




. sz L 






— <r- 





O (X 



o 



2 



zl-i 



I 



-O 



The change from \1/ to ^ removes the proportionality factor, which leads to (3.14). 



B.3. Proof of eigenvector property. We provide here a graphical proof of the eigenvector 
property of section 4.2. What follows is not rigorous because the VO construction becomes 
divergent when \q\ = 1; however, the proper proof, given in the text, is along the exact same 
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lines as this graphical proof, except with diverging prefactors removed. We shall therefore 
ignore prefactors in what follows. 

Recall the dual vertex operator (3.6). The properties (3.7) it satisfies can be depicted as: 



(B.3a) 
(B.3b) 



$ h (z)$* v (z) = S btV g 

i 



z . 


z , 








-< — 



9 



b=0 



Now use the first relation to rewrite the action of the transfer matrix T(z) on ^(z) in a 
similar way as in (4.5): 



o 



f \ 



oc 



o 



o 



Applying the qKZ equation (B.2) at i = L and ignoring the spectral parameter shift since 
s — 1, we find: 



Zl 



o 



ZL 



o 



oc o 



zl 

-( 



- 



(X o 



zl 
— *■ 



-O 



which is the desired eigenvector property. 

Appendix C. Properties of the spin-loop mapping 

We discuss here in more detail the linear map from the loop space to the spin space 
that is used in section 4.3. 

C.l. Case 1=1. This case is well-known, and we only review it briefly. tp\ = if : C[Ci] — > 
(C 2 )® L is defined explicitly in 4.3.1. Both spaces come equipped with standard bases, indexed 
by Cl and {0, 1} L respectively, and its matrix is given by (4.6); its entries are powers of —q, 
and therefore ip is well-defined for any q £ C x . Up to an overall power of q, these are 
nothing but maximal parabolic Kazhdan-Lusztig polynomials, see e.g. [ '] for a discussion 
in a related context. 

There is a partial order on {0, 1} L : e < e' iff Yl]=i( e j ~ e 'j) — f° r alH = 1, . . . , L. Denote 
e = {1, 0, . . . , 1, 0} and Cl = {e G {0, 1} L : e > e } (Dyck words). There is a natural 
bijection between C L and Cl which to n E Cl associates e such that 



1 if 7r(i) > i 
if < i 



i = l,...,L 
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The claim is that the submatrix of where indices are restricted to Cl x C l , is upper 
triangular w.r.t. the partial order above (modulo the bijection). This is fairly obvious in 
the graphical interpretation where one orients the arcs: the diagonal elements of the ma- 
trix correspond to all arcs pointing right, and each subsequent change of orientation of an 
arc clearly increases the partial sums X]j<i e «- Furthermore, the diagonal elements are all 
(~q) L ^ 2 - Therefore, the matrix has maximal rank for all q E C x , and cp is an isomorphism 
of C[£l] onto its image. 

For L = 2, the unique link pattern is sent onto the state (— q)~ 1 / 2 v w + (—qY^Voi, which 
is checked explicitly to be Lq-invariant. One can then proceed inductively by noting that 
any link pattern possesses an arc connecting neighbors, applying the L = 2 observation to 
it and then removing it. Therefore all link patterns are sent into the Lq-invariant subspace 

{{c 2 r L ) mv of (c 2 r L . 

For generic q, the dimension of the invariant subspace is the same as for sl(2) and is 
well-known to be the Catalan number L!/((L/2)!(L/2 + 1)!) (L even), which enumerates 
Dyck paths or link patterns. Therefore, for generic q, ip is an isomorphism of C[£l] onto 

((C 2 )^) 

inv ■ 

C.2. Fused case. Let us first define the projection map p = : (C 2 )® e — > V == C i+1 . We 
use induction on £: p^ = 1 and 

v (k+i) = (k) ( - [ , 1 k - r k \ ik) 

/0 0\ * 

where e fc = ( T -q- 1 o J ' m ^ ne ^°°P l an g ua g e , e k is simply * — * . From its defini- 
Vo o o o/ k,k+i i ^ 

k k+i 

tion, it commutes with the Lq-action, and in fact the claim (see e.g. [23]) is that it is exactly 
the projection onto the maximal spin £/2 irreducible subrepresentation. 

Now define (fe from (p by projection, as in (4.7). It is well-defined when the projection 
is, that is, when q 2r ^ 1, r = 2, ...,£. By definition, the (pi(ir) live in the [/i-invariant 
subspace {V® L )i nv of V® L . Furthermore, borrowing an argument from [32], we note that 
linear independence of the ^>e(^) is equivalent to linear independence of the ip(n), because 
ipg{ir) = V 9 ( 7r ) + J2 P ec iL -Ct L c 7r,p^(p)- But the latter was proved in the previous section. 

It is not hard to show once again that dim(V® L )i nv = #Ce,L for generic q (a short proof 
is to use a bijection - generalizing the one of the previous section - between fused link 
patterns and Bratteli diagrams for spin £/2, the increments of the paths being now given as 
the difference of numbers of opening and closing arcs in a given group); equivalently, using 
the crystal limit q — > one can conclude directly that the map above is an isomorphism 
between C[£ £ , L ] and {V® L ) inv for generic q. We conjecture that it remains an isomorphism 
for q = —e~ t7T ^ e+2 \ though we do not have a rigorous proof of that. 

C.3. Wheel condition. We prove here the wheel condition of Theorem 3.4. 

We work inside (C 2 )®^. We fix a subset of consecutive 2£ indices, say £(i — 1) + 1, . . . , £{% + 
1), and consider the projector Pj&i+i (j = 0, . . . , £) onto the spin j subrepresentation of the 
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action of U q ($l(2)) acting on these 2£ spaces. Also, given a link pattern 7r G Cl,i- denote 
2cj ) j + i(7r) the number of connections between the sites £(i — 1) + 1, . . . ,£(i + 1) and the 
outside. 

Then we have the following lemma: j > Cj ) j + i(7r) implies Pj. )iji+ iipi(TT) = (and the 
Pj-i t i + iifi(Tr), where it runs over link patterns with c^+i = j, are linearly independent). 
The equality is easy proved by inductively removing the "little arches" connecting sites 
£(i-l) + l,...,£(i + l). Explicitly, suppose n(k) — k+1, £(i — i) + l <k< k + 1 < £(i + l), 
and consider the vector space 

V = (C%_ 1)+1 ® ■ ■ • ® (C 2 )^ ® {{-q)- 1/2 |tl) - (-g) 1/2 |4t»*,W.i 

® (C 2 ) fc+2 ® • - • ® (C 2 )^ +1) 

Note that ^(tt) G <8) K g> W right , where W left S (C 2 )®^" 1 ) and W rig ht = (C 2 )®^-*- 1 ) 
correspond to the sites outside [£(i — 1) + l,£(i + 1)]. Since the vector at sites k, k + 1 is 
a singlet of C/ g (sl(2)), V is isomorphic to (C 2 )® 2 ^ -1 ^ as a representation space of Z7 9 (s[(2)). 
Iterating the process, we find that G W\ e ft ® V <g> Weight as above, where V is a subspace 

of (C 2 )® 2£ which is isomorphic to (C^^m+iO). But we know explicitly the decomposition of 
the latter into irreducible representations, and in particular that the highest possible spin is 
Cj i j + i(7r), hence the equality of the lemma. For the linear independence, consider the unique 
spin state obtained from ir which has all outside connections of the form i.e., by orienting 
all lines to the outside "outwards". 

Next, we note that since the projection n^=iP« cannot increase the number of connections 
to the outside, then the same statement can be made about the "fused" vector <pe(ir), n G Cl,i 
(noting that in the case of fused link patterns, the connections to the outside are shared 
equally between groups i and i + So, 

Pj-i,i+i<Pi(ir) = 7T G £ L; £, j > c iti+1 (n) 

(and the Pj-^i+iip^Tr), where n runs over fused link patterns with c^+i = j, are linearly 
independent). 

From the coproduct of U q (si(2)) it not hard to see that Pj.^+i commutes with its action 
on the whole of (C £+1 )® L , and in particular leaves invariant the image of ifi£. This implies 
that Pj.i t i + i(pt(Tr) is a linear combination of <^(7r'). 

Now writing P J j.j ) j + iP 3 -.^ +1 <^(7r) = for j < £, and expanding Pj-i^xfti^) in (pi(n f ) as 
above, we find that no link patterns tt' with 0^+1(71"') = £ can appear in the expansion. 
Repeating the argument with P e _ 1 . iti+1 P j . iji+1 (p l (Tr) = 0, etc, until P j+ x-i,i+xPj;i,i+i(pe(^) = 0, 
we conclude that 

ImP^j+i o(p e C ((fe(ir) : Cu +1 (n) < j) 

Now consider ^ := \&(°)(. . . , z», z i+ i, z i+2 , . . .) (the case k > will be discussed separately) 
with Zi = q 2( ~ a+b 'z, z i+ \ = q 2b z, z i+2 = z, a, b > 0, a + b < £ + 2. Plugging these equalities into 
(3.13) and using the form (3.12) of the F-matrix, we conclude that ^ is a linear combination 
of Pj-i t i + i^ , j < a, and a linear combination of Pj-i + x ) i+2 i & , j < b. 

As mentioned in sect. 4.3, is [^-invariant. We can therefore apply the reasoning above 
and conclude that 

* G ((M 71 ") : Q,i+i(7r) < a, Q +M+2 (7r) < b) 
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For any such ir the group of sites i + 1 has at least £ + 1 — a connections to the sites i, and 
£ + 1 — b connections to the sites % + 2; but £ + 1 — a + £ + 1 — b > £ which is contradictory. 
Therefore ^ = 0. 

This is the special case i\ — i, i<i — i\ + 1, is — i\ + 2 of (3.16). In order to obtain 
the general case, one simply needs to apply the exchange relation (3.13) repeatedly, to move 
(i,i+l, i+2) — > (ii, i 2 , is); assuming that all other arguments of are generic, all .R-matrices 
involved in this process have nonzero denominator and are therefore well-defined. We have 
thus proved Theorem 3.4. 

The case of vj/w, k > 0, can be treated similarly by noting that the reasoning above 
is purely local and cannot depend on the twist. Explicitly, one can build a [/i-invariant 
element of C fc+1 <8> (C 2 )® eL <g> C fc+1 given by (x\Q(zi) . . . $(^2n) \v) where (x\ and \y) run over 
the [/i-representation generated by (k\ and \k), each isomorphic to C k+1 . It satisfies the 
wheel condition, and therefore \?' fc ) does. 

Remark: the connection between the projectors Pj and link patterns can be made more 
explicit by introducing generalized Temperley-Lieb operators Ej (denoted e"~i' in [32]), 
which are best described graphically as 



Ej = P1P2 




Note that the Ej are (up to normalization) a family of non-orthogonal projectors (as opposed 
to the Pj which are orthogonal to each other). 

Comparing the expression of the i?-matrix (3.12) with the following expression (correcting 
a sign mistake in Eq. (2.8) of [32]): 

we see that there is a triangular change of basis between the Pj and the Ey. Pj = J2j><j a j,j'Ej' 
with ctjj 7^ 0. 

Appendix D. Polynomiality 

In this appendix we prove that the expression (3.10) is a polynomial. As a bonus we 
get that the coefficients are Laurent polynomials in q up to some known factor. Recall the 
expression: 

2n I „ in 



^)-7&'n^ n n(^--.)f--/nS 

i=l Ki<j<2nr=l J J i=l 

n 



(D 1) 1=1 X - i<J 

w . PyW(Wl,...,W tn ) 

X ■ ■ 



. Wj - q 2 Wi Ylj<eM eWi - z i) Uj>eS Wi ~ ^ z i) 
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where the contour integration is made in the following way: we start by integrating wg n 
arround the points q e ~ 2k z s , for all 1 < % < 2n and < k < £; we proceed by doing the 
same with ru^-i, W£ n _ 2 , up to w\. By Cauchy's residue theorem this is the same 
as replacing {wi, . . . , wi n } — > {q e ~ 2kl z si , . . . , q e ~ 2ken z S(n }, removing the correspondant poles 
(wi — q e ~ 2kt z s .)~ 1 and summing over all possible combinations. Notice that, due to the 
presense of the Vandermonde polynomial we can not repeat a pole. Also, q e ~ 2k z s is only 
possible, for k > 1 if q l ~ 2k+2 z s appears to the right. 

We want to prove that there is neither any pole of the kind (z s — q a z T )~ 1 nor of the kind 
z~ l . In order to prove it we will look in the details of the computation. We can ignore 
all integrations except the ones that concern z r and z s , with r < s. There is a subset of 
the integration variables {wj 1 , . . . ,Wj m } which is replaced by q e ~ 2h z r and q e ~ 2k z s , the others 
being replaced by other variables z t . 

Fix m. We will represent it diagramatically, representing all poles in z s by a blue circle 
and the ones in z r by a red square. For example, we represent 

{w 1 , w 2 , w 3 , Wi} {q e ~ 4 z s , q e ~ 2 z s , q e z r , q e z s } by • • ■ • 

where we simplified the indices of the integration variables. 

Let ipd,e be the relevant part of the integral, where d stands for an arbitrary diagram with 
m circles and squares, 



in 



dWi -p-r Wi — Wj T-r 1 



^dA Z r,Z s ) = f ... fY[—^Y\ -2 II V ' W ' \ 

J J -M; 27T« t^t Wi - q 2 Wj - L .- L [Wi - q E z r ){Wi - q l z s ) 
(D.2) l ° 1 t 



2., 



where the integration is performed according to the diagram d. We restrict e to the relevant 
indices: {e h , . . . e jm }. 

The following statement is necessary to the proof: 
Proposition D.l. The Macdonald polynomial P v (k) satisfies: 

P Y(k) (z, q 2 z, . . . , q 2 i- 2 z, w j+1 , . . .) = z ™ x ^ k ^ x (. . .) 

where j < £. 

Recall that Py) is obtained from Wi) (ki,...,kt \ e^^J 7 \ki, . . . , ke), which can be seen 
as a sum over all parafermionic decompositions. Chosing j variables {z, q 2 z, . . . , q 2 i~ 2 z\ 
forces the parafermions to be in distinct modes, according to (A. 23). Each parafermion 
f l (z) generates a factor z~ ki . The result follows from minimizing the power of z on the 
expression. 

We split formula (D.2) into two contributions, the first row and the second row. All the 
poles of the second row are outside of the contour integral. 
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Proposition D.2. For any diagram d, the second row of formula (D.2) is a polynomial in 
z r and z s . Moreover, if multiplied by j, ^ , where m r (and m s ) is the number of Cj 
egna/ to r (respectively s), it is a Laurent polynomial in q. 

Notice that m r < b r (and m s < b s ), the equality being only possible if we chose all W\ 
such that €i = r (or s) to belong to our subset of variables. 

Proof. By choice r < s. Pick a random diagram, for example 




r s 



where the vertical lines split the variables W{ in five regions q < r, q = r, r < < s, = s 
and ej > s. The product (ifj — q e z r )(wi — q e z s ) vanishes whenever there is a blue circle on 
the region q > s; and it vanishes also if there is a red square on the regions q > r. Then 
this example will vanish. 

We pick now a non vanishing example: 




r s 



The term (g £ u?j — z s )~ l only appears in the region = s, exacly m s times, when we replace 
Wi by the sucessive q e ~ 2p z s we get (q 2e — . . . (q 2£ - 2m s+ 2 _ l)- 1 ^"" 18 . There is at least 
m s variables z s in P v (k) and we can apply proposition D.l, therefore the power in z s is at 

least z^zf a ' x ^ rna ~ k ' ^ z~ ms = <2 ™ ax { ' fc ~ m 'd. ^he term (g £ w, — z r )~ l appears in all regions such 
€i > r. Let m' r be the number of red squares in the region e« = r. Then we can do the same 
analysis for the squares. The circles will generate a term (q 2e z s — z r ) . . . (q 2e - 2m +2 z s — z r ), 
where m! is the number of circles in the regions e.; > r. This term cancels with the product 

Il P (<l 2P Zs-Zr). □ 

In general, the first row of the formula (D.2) will give a rational function. 

Proposition D.3. The sum over all possible diagrams with m nodes of i/Jd,e( 
polynomial in z r and z s . 

Proof. From proposition D.2 we know that 

dWi -p-r Wj - Wi y-r Qd(z r , Z s ) 



i 2ttz jj: Wj - q 2 Wi A . A (w { - q e z r )(wi - q £ z s ) 

where Qd{z r , z s ) is a polynomial depending on the diagram d. We ignore the dependence in 
e. For convenience we allow diagrams with more than I circles (or squares), in which case 
we set Qd to zero. This is naturally done by using the Macdonald Polynomial Py t ■ 

On the other hand, Qd( Zy^ ZgJ Cclll be seen as the evaluation of a function, i.e. Qd(z r , z s ) = 
Q(z r , z s ; z^\ . . . , z^ 71 '), where is given by the diagram d. We will omit the explicit 
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dependence on z r and z s , and for the analysis we will consider that the are independent 
variables. Then, this is a rational function of the form 

R(z^ z^) 

(D.3) Q(z« . . . , z (m) ) = 1 j 



where the polynomial R and Q. depend on z r and z s . 

The diagram • — ■ — • — • corresponds to the following expression: 

, / \ Q{q i ~ 4 Zs,q £ z r ,q e ~ 2 Zs,q e z s ) 

l/Jd{Zr,Zs 



(q^Zs - q e - 2 z r )(q £ z r - q e ~ 4 z s )(q e ~ 2 z s - q e z r )(q e z s - q l z r ) 

The denominator can be constructed by a simple rule: we read the diagram d from right to 
left, at each node we can put either a circle or a square (corresponding to some q a z s or q b z r ); 
if we choose a circle (resp. a square) we add the pole (q a z s — q b z r )~ l (resp. (q b z r — q a z s )~ 1 ). 

Now we sum over all 2 m possible diagrams. Let d\ and d 2 be two diagrams that only differ 
in the first element, i.e. d\ = md^ and d 2 = •d^\ where d^ is some diagram with m — 1 
nodes. Then: 

1 Q dm ~ Qmdm 



i>m d W (Zr, Z s ) + 1p. d (i) (z r , Z s ) 



q a z r - q b z s 



where the three dots are replacing the commun m — 1 poles. Using the formulation (D.3) 
and setting q a z r = z^ and q b z s = z'^\ the numerator is given by 

R(z^\...) i?(z'«,...) 



ru n jt {*® - oj \u n (z^ - n h (z' il) - 

which vanishes if z^ 1 ' = z'^\ Then the division by (z^ — z'^) is well defined and we can 
define the polynomial: 

q(1) . = QmdW - Q»dW 

y As q A r 

where the open circle indicates that we already performed the sum. This can be written in 
the rational function perpesctive: 

Q {1) (z {1) , z' {1) ; z {2) . . . ,z {m) ) - 



Notice that this is a symmetric function on z^ 1 ' and z'^\ 

Now we sum over the second node, let d\ = omd^ and d 2 = °»d^ 2 \ where the open circle 
indicates the first node where we already performed the sum. Then: 



1 Q V> _qW 

Wo. d m\Zr,Z s ) +iPomdi*){Z r ,Z s ) = ? T - ^ 



q a z s - q°z 
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We consider now that z^ 1 ' and z'^> depend on the two variables and z'^ 2 \ then the 
numerator is given by: 

i?w( ? -y 2 ), 2 '( 2 ^( 2 ),...) 



i?(l)( z (2) ) 5-2^(2)^/(2)^ _ 



n^ (2) - c*)(?-w - a) n> /(2) - 

Once again this vanishes if z^ = z ,<y2 \ then the division by (z^ — z'^ 2 ') is well defined and 
we can define a new polynomial: 

or in the rational function perspective: 

w 1 • • ' n^ (2| -4)(^ 2 »-a)n,> 2 a,(^»-c„) 

which is once again symmetric in z^ and z'^ 2 \ 
Repeating the process m times we obtain: 

ip(z r ,z s ) = Q {m) (q e z r ,q £ z s ) 
which is a polynomial. □ 

This concludes the proof. Note that if Qd is a Laurent polynomial in q for all d, then 
ipd will be also a Laurent polynomial in q. Then: 

Corollary D.4. The coefficients of ([^]!)~ n IU^]^^... b^ n { Zl i ' ' ' > Z2n ^ are Laurent polyno- 
mials on q. 
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